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As a secondary payload poised atop a Delta rocket, the Danish geomagnetic research 
satellite Oersted was launched February 23, 1999 from Vandenberg AFB in California. 
The principle aim of the Oersted mission is to accurately map Earth’s global magnetic 
field, a field which has changed significantly in the twenty years since NASA’s 
enormously successful Magsat mission. After attaining a near polar orbit with altitudes 
ranging from 638 to 849 km, the mini-sat class Oersted spacecraft extended an 8 meter 
long gossamer mast carrying both vector and scalar magnetometers. These instruments 
measure both the direction and strength of the ambient magnetic field. The orientation of 
the measured field relative to the fixed stars is determined using an effectively non- 
magnetic Star Imager bonded to the vector magnetometer. The Oersted Initial Field 
Model is derived from both scalar and calibrated, oriented vector magnetic data using the 
time-tested technique of spherical harmonic analysis. The thousands of measured data 
fitted by this main field model were carefully selected to reduce contributions from 
electric currents flowing in Earth’s partially ionized upper atmosphere and from magnetic 
storms caused by the fluctuating solar wind. Moreover, the weighted least squares fit to 
these select data reflects the anisotropic residual covariance needed to account for small 
errors in the observations, notably in the angle of rotation about the bore-sight of the 
single Star Imager. As expected, the model fits the measured field to about one part in 
ten thousand. Although the Oersted Initial Field Model includes Gauss coefficients 
through degree and order 19, we only recommend use of coefficients from this initial 
field model up to degree 14 at most. 
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Abstract. Magnetic measurements taken by the 0rsted satellite during geomag- 
netic quiet conditions around January 1, 2000 have been used to derive a spherical 
harmonic model of the Earth’s magnetic field for epoch 2000.0. The maximum 
degree and order of the model is 19 for internal, and 2 for external, source fields; 
however, coefficients above degree 14 may not be robust. Such detailed models 
exist for only one previous epoch, 1980. Achieved rms misfit is 2 nT for the scalar 
intensity and 4 nT for the vector components perpendicular to the magnetic field. 
This model is of higher detail than the IGRF 2000, which for scientific purposes 
related to the 0rsted mission it supersedes. 


Introduction 

Twenty years after the Magsat mission, the 0rsted satel- 
lite was launched on February 23, 1999 in a near polar or- 
bit with an inclination of 96.5°, a perigee at 638 km and 
an apogee at 849 km. The principal aim of the 0rsted mis- 
sion [Neubert et al , in preparation] is to accurately map the 
Earth’s magnetic field. 

The goal of this paper is to present an accurate “snapshot” 
of the geomagnetic field at epoch 2000.0 based on geomag- 
netic quiet days around January 1, 2000. To reduce contam- 
ination of the lower degree expansion coefficients by spatial 
aliasing, the analysis was performed through degree and or- 
der 19 for the internal field and 2 for the external field. How- 
ever, we only recommend use of coefficients up to degree 14 
at most. 

0rsted data from May to September 1999 were used to 
derive the IGRF 2000 [Olsen and Sabaka , 2000]. How- 
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ever, the satellite was still in the commissioning phase before 
September 1999, and the measurements used for the present 
study are more accurate. 

Data selection and pre-processing 

Data from geomagnetic quiet conditions between Decem- 
ber 18, 1999 and January 21, 2000 were selected according 
to the following criteria: Kp < 1+ for the time of obser- 
vation, Kp < 2o for the previous three hour interval, |Dsf| 
< 10 nT and \d(Dst)/dt\ < 3 nT/hr. To reduce contribu- 
tions from ionospheric currents at middle and low latitudes, 
only night-side data (local time about 22:00) were selected. 
Vector data were used for dipole latitudes |A d i p | < 50°, and 
scalar data for |A d i p | > 50° or if attitude data were not 
available. To further reduce contributions from polar cap 
ionospheric currents, only data for which the dawn-dusk in- 
terplanetary magnetic field component was weak (|B y | < 
3 nT) were used. An equal area distribution was approxi- 
mated by decimating the data to measurement times at least 
30s/ sin 6 apart, where 0 is geocentric co-latitude. All orbits 
were visually inspected and those suspected of contamina- 
tion by external currents were removed. Attitude informa- 
tion over the South Atlantic Anomaly (SAA) was sparse due 
to radiation effects on the star imager, resulting in a paucity 
of vector data there. Since this region is close to the dip- 
equator (the thin line in Figure 1) where vector data are 
mandatory to avoid the Backus effect [Stem and Bredekamp , 
1975], 66 data points during four orbits on the geomagnetic 
quiet days October 9 and November 26, 1999 were added to 
constrain field direction in this region. Figure 1 shows the 
distribution of the 2148 scalar data points and 3957 vector 
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Figure 1. Data distribution. Scalar measurements are shown 
by small and vector measurements by larger symbols. Open 
circles present additional vector data to fill the gap in the 
SAA. 

triplets used for the model. 

The attitude accuracy of the 0rsted star imager (SIM) is 
anisotropic: determination of the SIM bore-sight direction 
(“pointing”) is more accurate than determination of the ro- 
tation around the bore-sight (“rotation”). This results in rel- 
atively more noise in the rotation angle, which is also more 
sensitive to distorting effects like instrument blinding (for 
instance by the moon). This attitude anisotropy results in 
correlated errors between the magnetic components, which 
should be taken into account when deriving field models 
[Holme and Bloxham, 1996; Holme, 2000]. 

Let n be the unit vector of the SIM bore-sight, and let 
B be the observed magnetic field vector. B and n define 
a new coordinate system (provided they are not parallel), 
and the magnetic residual vector £B = (SB b ,6B ± ,8B 3 ) 
can be transformed such that the first component, SB b , is in 
the direction of B, the second component, 8B±, is aligned 
with (nxB), and the third component, 6B 3 , is aligned with 
B x (nxB). The last two components are perpendicular to 
the magnetic field. In this coordinate system, the errors on 
the different field components are uncorrelated. 

Let xj) be the attitude error of the bore-sight (pointing er- 
ror), x that around bore-sight (rotation error) and let a be 
the (attitude independent) error of the scalar intensity. Con- 
sidering only linear terms in xjj and the component SB b 
is independent of attitude errors. 5B 3 is affected only by 
pointing errors, whereas SB± is influenced by both pointing 
errors and rotation errors [Holme, 2000], Since the rotation 
uncertainty is believed to be the main error source, the com- 
ponent Bx is more contaminated than the two other compo- 
nents and should be down-weighted. 

As an example, Figure 2 shows the residuals of the 
nightside part of orbit # 4343. The upper panels present 
residuals ( SB r , 6B S , 5B $) as a function of colatitude 6\ the 
lower panels present the residuals as (SB b ,8Bx,5B 3 ) . The 
noise is clearly spread over all three components in the 



Figure 2. Residuals (observations minus values predicted by 
the model of Table 2) as a function of co-latitude 6 for the 
nightside part of orbit # 4343, December 22, 1999, 01:21 - 
02: 1 1 UT. Kp = Oo and Dst = +9 nT. 

0 5B r ,8Bo,8B 0) system, but is concentrated in 8B± in the 
( 5B b ,5Bx,8B 3 ) system. 8B 3 remains slightly more noisy 
than 5B b due to pointing uncertainty and field-aligned cur- 
rents. In particular, 8B 3 signatures at 6 « 25° and 170° 
are caused by auroral field-aligned currents. The gap in the 
vector components between 6 = 100° and 135° is due to 
moon-blinding of the star imager. 

Most all data selected span 35 days during which the field 
changes by up to 20 nT due to secular variation (SV). It has 
been decided to account for this change by propagating all 
observations to epoch 2000.0 since this a) reduces the model 
misfit of SB b and 5B 3 by 10% and 6%, respectively (that 
of 8Bx is not changed); and b) reduces the magnetic power 
Rn for all degrees above n = 12 (mean reduction: 17%). 
The latter indicates that SV between neighboring measure- 
ments taken at different epochs will, if not corrected for, pro- 
duce spurious high degree signals. Two different SV models 
were used: an updated version of the IGRF 2000 SV model 
[Macmillan and Quinn, 2000] and a model generated from 
0rsted data spanning over 9 months. Negligible difference 
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was found, so it was decided to use the 0rsted SV model for 
internal consistency. This model is still in development, and 
will be the subject of a future publication. 

Model parameterization and estimation 


The magnetic field B = -grad V can be derived from a 
magnetic scalar potential V which is expanded into terms of 
spherical harmonics: 
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a = 6371.2 km is the mean radius of the Earth, (r, 0,^) 
are geocentric coordinates, P™ (cos0) are the associated 
Schmidt-normalized Legendre functions and [g™, h™) and 
(g£\$™) are the Gauss coefficients describing internal and 
external source fields, respectively. The coefficients 
and s\ account for the Ds£-dependent part of the external 
dipole. Its internal, induced, counterpart is represented via 
the factor Qi = 0.27, a value found from Magsat data by 
Langel and Estes [1985]. The model has 410 parameters 
(399 static internal, 8 static external, and 3 coefficients of 
Dsf-dependency). 

The coefficients are estimated by an iterative least- 
squares fit, minimizing e r C^ 1 e where the residuals e = 
dobs “ d mo< j are the differences between observations and 
values predicted by the model, and C d is the a-priori covari- 
ance of residuals due to data errors and fields left unmodeled. 
To account for the anisotropic attitude accuracy, C d was as- 
signed the form of Eq. 18 of Holme and Bloxham [1996], 
with a — 2.25 nT, = 10" and x = 60" (these values are 
justified by the a-posteriori model misfit, pre-flight instru- 
ment error estimates are 3" and 20" for SIM and 0.2 nT for 
the vector magnetometer). 

When solving the least-squares problem, three iterations 
were found to be sufficient. Outliers were removed after the 
second iteration; as outlier selection criterion we have used 
15 nT for SBb and SB$ but 30 nT for 5B±. If one of the 
component residuals was above its threshold value, all three 
components were removed. 

Results and Discussion 

Number of data points fitted, residual means and rms mis- 
fits of the model are given in Table 1. The anisotropy of rms 
misfit in ( SB b , 5Bj_ , SB 3 ) components is much larger than 
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0.04 
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8.38 

6B 3 
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2.62 

6B r 

3957 

49 

0.63 

4.79 

SBg 

3957 

49 

-0.20 

5.58 

SB# 

3957 

49 

-0.36 

5.21 


Table 1 . Number N to t of data points, number A^ out of re- 
moved outliers, means, and rms misfits (in nT) for the dif- 
ferent components. "Polar” denotes data with |A dip | > 50°. 

in spherical components. Figure 3 shows the data residu- 
als as a function of dipole co-latitude. The largest residuals 
in SB b are found in the southern polar cap (0 dip > 170°) 
and are attributed to ionospheric currents in the summer po- 
lar cap. Electrical conductivity is smaller in the northern 
(winter) polar cap ionosphere due to absence of solar irradi- 
ation, and therefore the SBb residuals in the northern polar 
cap (0 dip < 10°) exhibit a smaller scatter than their counter- 
parts in the southern polar cap. However, contributions from 
ionospheric polar electrojets are present even in the northern 
auroral zone (0 dip « 23°). 

The large rms and mean values for 6B± (present even 
when the data are weighted isotropically) confirm the neces- 
sity to downweight this component using the covariance ma- 
trix C d . A resolution analysis shows that B± resolves only 
6% of the model parameters if the data are treated in this 
way, but 24% if data errors are assumed to be isotropic. The 
corresponding values for B b (B 3 ) are 38% (27%) and 27% 
(23%), whereas the resolving power of F is 30% and 27%, 
respectively. The model fits the data with a (normalized) chi- 
squared misfit of 1 .01 which is close to unity, consistent with 
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Figure 3. Residuals as a function of dipole co-latitude. 
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the weighting (a-priori data errors <7, rp, x) being correct. 
The model coefficients are given in Table 2. Only internal 
coefficients up to n = 14 are listed; the complete model is 
available at www.dsri.dk/OerstedyField_models/OIFM/. 

Experiments with various truncation levels of the spher- 
ical harmonic expansion gave the largest changes in the 
southern polar cap (0 dip > 175°), centered around the mag- 
netic pole. This indicates that contributions from iono- 
spheric current systems in the summer hemisphere are prob- 
ably present in the data and model in spite of our attempt to 
minimize external current contributions by careful data se- 
lection according to geomagnetic indices and IMF B y . 

To address these contributions will likely require inclu- 
sion of more satellite data, ground data, and perhaps co- 
estimation of ionospheric field parameters. Doing so is be- 
yond the scope of our model. 

The model presented here provides a firm basis for stud- 
ies of the ionospheric, magnetospheric, lithospheric and core 
fields. The initial model from 0rsted, the first magnetic 
mapping mission of the “International Decade of Geopo- 
tential Research”, will also aid interpretation of the addi- 
tional, continuous, high precision measurements of Earth’s 
time-varying geomagnetic and gravitational fields acquired 
by forthcoming missions like Champ and S AC-C. 
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Table 2. Expansion coefficients of internal (g™, h™) and 
external ( q ™ , s™) contributions, in nT. q™ and s™ present 
the Ds£-dependent part of the external coefficients. 


